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Abstract 

The quantum search problem is an important problem due to the fact 
that a general NP problem can be solved by an unsorted quantum search 
algorithm. Here it has been shown that the quantum search problem could 
be solved in polynomial time on an NMR quantum computer. The NMR en- 
semble quantum computation is based on the quantum mechanical unitary 
dynamics that both a closed quantum system and its ensemble obey the same 
quantum mechanical unitary dynamics instead of on the pseudopure state or 
the effective pure state of the classical NMR quantum computation. Based 
on the new principle the conventional NMR multiple-quantum spectroscopy 
has been developed to solve the search problem. The solution information of 
the search problem is first loaded on the unitary evolution propagator which 
is constructed with the oracle unitary operation and oracle-independent uni- 
tary operations and used to excite the multiple-quantum coherence in a spin 
ensemble. Then the multiple-quantum spectroscopy is used to extract ex- 
perimentally the solution information. It has been discussed how to enhance 
the output NMR signal of the quantum search NMR multiple-quantum pulse 
sequence and some approaches to enhancing the NMR signal are also pro- 
posed. The present work could be helpful for the conventional high magnetic 
field NMR machines to solve efficiently the quantum search problem. 



1. Introduction 
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An efficient quantum search algorithm may be tremendously useful be- 
cause of the wide range of its practical applications, and especially due to 
the fact that a general NP problem can be solved by the unsorted quantum 
search algorithm [1]. It is believed extensively that a general NP problem 
can not be solved efficiently on a classical digital computer in polynomial 
time. But if there would be a polynomial-time quantum search algorithm, 
then a general NP problem would be efficiently solved by the algorithm in 
polynomial time on a quantum computer. In the past years a great effort has 
been devoted to searching for an efficient quantum search algorithm. On the 
basis of the quantum parallelism principle [2] Grover first proposed a quan- 
tum search algorithm with quadratic speedup over the classical counterparts 
[3, 4]. It has been shown that this algorithm is optimal [5, 6]. A large number 
of works of the Grover algorithm and its applications have been published in 
the past several years. However, the quantum search algorithms [3, 4] which 
are based on the pure quantum states are not polynomial-time algorithms 
and hence have not been proved certainly that they can move the frontier 
between solvability and intractability of complexity problems [6]. On the 
other hand, the experimental realization of quantum computation and the 
physical architecture of quantum computers have achieved a great progress 
in the recent years, although it is the long term to build up a commercial 
quantum computer. In particular, the experimental ensemble quantum com- 
putation based on the nuclear magnetic resonance (NMR) techniques [7, 8] 
has been extensively studied [9, 10] and a large number of the related works 
have been reported [11-14]. But the classical NMR quantum computation 
[9-15] working on the concept of the pseudopure state [9] or the effective 
pure state [10] suffers from many trouble problems [15, 16]. For example, 
people suspect whether or not the classical NMR quantum computation is 
real quantum computation and is as powerful as the quantum computation 
based on pure quantum states [16]. Since the effective pure states or the 
pseudopure states are isomorphic to the pure quantum states [9, 10, 12] the 
classical NMR quantum computation is at most as powerful as the original 
one based on pure quantum states and thus, can not be able to move the 
frontier between solvability and intractability. It has to give up many impor- 
tant characteristic features of a spin ensemble in order that with the help of 
the concept of the effective pure state or the pseudopure state the pure-state 
the pure-state quantum computation can work in the spin ensemble as well. 
For example, highly mixed states of a spin ensemble are not suitable for the 
classical NMR quantum computation. It is known that the quantum states 
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of the spin system in a spin ensemble are very highly mixed. Therefore, the 
classical NMR quantum computation can not give a satisfied answer whether 
or not the highly mixed states in a spin ensemble could be used in quantum 
computation and about the power question of quantum computation with 
the highly mixed states. 

Amazingly, recently it has been shown that an NMR quantum computer 
operating on the highly mixed states of a spin quantum ensemble can be 
exploited to solve efficiently the unsorted search problem [17] and the parity 
problem [18] in polynomial time, while these problems have been shown to 
be hard NP problems and can not be solved efficiently in classical compu- 
tation and even in pure-state quantum computation [19, 20]. These results 
show powerfully that (a) a general NP problem could be solved in polyno- 
mial time on a quantum computer; (b) quantum entanglement of a quantum 
system may not be the sole origin of quantum computation power and the 
quantum mechanical unitary dynamics may play an important role in solving 
these NP problems in polynomial time; (c) quantum computers can outper- 
form classical computers and ensemble quantum computation can be at least 
as powerful as the standard quantum computation based on pure quantum 
states. These quantum algorithms or sequences that operate on highly mixed 
states in a spin quantum ensemble could achieve an exponential speedup over 
the classical search computation and therefore they really move the frontier 
between the solvability and intractability of complex problems. Although 
in theory the NMR ensemble quantum computation on highly mixed states 
is successful to solve the unsorted quantum search problem and the parity 
problem in polynomial time, it still suffers from some problems and limits in 
practice. The most severe problem is that the output NMR signal intensity of 
these sequences reduces exponentially as the qubit number of the spin system 
in the spin ensemble [17, 18], making the NMR ensemble quantum compu- 
tation available only in spin systems with a small number of qubits. There 
may be two schemes to extend the scale of the NMR quantum computation 
so that it could solve the search problem with a larger qubit number. One 
is to enhance the initial spin polarization of a spin ensemble by a variety of 
possible spin-polarization enhancement techniques [7, 8, 21, 22, 23]. Another 
is to design new quantum search algorithms or sequences which consist of a 
polynomial number of the oracle unitary operations and oracle-independent 
unitary operations since the output NMR signal could be increased by many 
calls of the oracle unitary operations in the sequences. In this paper ac- 
cording to the second scheme new quantum search sequences on the NMR 
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ensemble quantum computers are proposed to solve the search problem with 
the help of the multiple-quantum operator algebra space formalism [24, 25], 
and particularly the NMR multiple-quantum spectroscopy [7, 8] has been 
developed to search for the unknown marked state, i.e., the solution of the 
search problem. 

There may be two possible schemes for quantum computation making a 
transition from a pure-state quantum system to its quantum ensemble. One 
is the classical NMR quantum computation. It is based on the concept of the 
effective pure state [10] or the pseudopure state [9] which is isomorphic to 
the pure quantum state. Then, the quantum computation on the pure state 
quantum systems allows to apply to their quantum ensembles with the help 
of the concept of the effective pure state. In the frame of the classical NMR 
quantum computation information of quantum computation is mainly car- 
ried by the effective pure states. Another is the new NMR ensemble quantum 
computation. It is based on the fact that both the closed pure-state quan- 
tum system and its quantum ensemble obey the same quantum mechanical 
unitary dynamics when decoherence effects in the quantum system and its 
ensemble are negligible [17, 18]. Then in the frame of the NMR ensemble 
quantum computation the information of quantum computation is usually 
loaded on the unitary evolution propagator that describes the quantum me- 
chanical unitary dynamics of the quantum system and its ensemble. The 
highly mixed states could be used in the ensemble quantum computation be- 
cause the new ensemble quantum computation allows the independence each 
other between the initial input state and the quantum circuit in a quantum 
algorithm or sequence. Therefore, the quantum search process based on the 
new ensemble quantum computation is that the information of the unknown 
marked state of the search problem is first loaded on the quantum mechanical 
unitary evolution propagator, then one uses various experiment and measure- 
ment techniques such as the NMR multiple-quantum spectroscopy to extract 
the solution information of the search problem. 



2. The simple quantum search sequences on a spin en- 
semble 

First of all, the explicit form of the oracle unitary operation of the quan- 
tum search problem is introduced below. In a quantum search algorithm the 
basic oracle unitary operation Uf is defined as [26] 
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U f \x)\a) = \x)\a@f(x)) (1) 
where the quantum state \x)\a) is an arbitrary computational basis of the 
whole quantum system which consists of the work qubits I and the auxiliary 
qubits S, the quantum states \x) and \a) are the computational bases of the 
work qubits / and the auxiliary qubits S, respectively, and the direction sum 
symbol denotes addition modulo 2. The function f(x) in the quantum 
search algorithm is defined over the range < x < N — 1 (N = 2 n ) by 
f(x) — 1 if x is a solution to the search problem (x = s), which corresponds 
to the marked quantum state \s) of the work qubits /, and f(x) — if x is 
not a solution to the search problem (x ^ s). 

When the basic oracle unitary operation Uf is applied to the quantum 
state | a;) 1 5) where the quantum state \S) of the auxiliary qubits S is chosen 
as the superposition: \S) = ^QO) — |1)) it follows from the definition (1) of 
the oracle unitary operation Uf that [26] 

U,\x)\S) = exp(-i7r/(x))|x)|5) = { = * . (2) 

This shows that only the marked quantum state \s) is inverted in phase 
or changes sign of its amplitude but any other state \x) keeps unchanged 
when the function values f(x) are simultaneously evaluated once over all x 
values, that is, the oracle unitary operation Uf is performed once on the 
quantum system. One can see from Eqs.(l) and (2) that the oracle unitary 
operation Uf is independent of any quantum state \x) of the work qubits / 
but related to the quantum state \S) of the auxiliary qubits S. It must be 
paid much attention to this point [17, 18] when a quantum search algorithm 
constructed with the oracle unitary operation Uf and oracle-independent 
unitary operations is extended to a quantum spin ensemble from its pure- 
state quantum spin system. For convenient design of a quantum search 
algorithm or sequence, an auxiliary oracle unitary operation U Q is introduced 
below. This oracle unitary operation is expressed as U Q {9) = Uj 1 Vs(0)U / 
[17] whose unitary operation on any quantum state of the quantum system 
is defined by 

U o {0)\x)\m) = ^V 5 (0)|x>|Oe/Or)>|l) 

= exp(-^(/(a;),l))|x>|0>|l> (3) 
where the conditional phase-shift unitary operation Vs(0) [27] acts only on 
the auxiliary quantum state \S) = \a)\b) and is defined by 

^)|a)|l) = exp(-^(a,l))|a)|l) = { ^ ( ~** >|1} ^ = J ■ (4) 
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Note that the basic oracle unitary operation Uf has the property: Uj = E 
(unity operation), as can be seen from Eq.(l). Then the oracle unitary 
operation U can also be expressed as U {9) = UfVs(6)Uf. 

The quantum parallelism principle [2] provides the possibility for quan- 
tum computation outperforming the classical computation. By the quantum 
parallelism all the function values f(x) can be evaluated simultaneously over 
all the input computational base \x) of the quantum system when the oracle 

N-l 

unitary operation U acts once on any superposition |^) = E a x \x)\S) of 

x=0 

the quantum system, 

N-l 

Uo(0)\*) =E exp(-i96(f(x),l))a x \x)\S) 

N-l 

= exp(-i9)a a \s)\S)+ E a x \x)\S) (5) 

x=0,x^s 

where the definition of the function f(x) is used: f(x) — 1 if x — s; f(x) = if 
x 7^ s. To analyze conveniently the unitary evolution process of the quantum 
system and its ensemble during the action of the oracle unitary operation U Q 
one needs to use the selective phase-shift operation C s (9) [17] to express 
explicitly the oracle unitary operation. The selective phase-shift operation 
C r {9) is a diagonal unitary operator which is only applied to the work qubits 
/. It can be expressed in the exponential form 

C r (9) = exp(-i9D r ) (6) 
where the diagonal operator D r is related to the quantum state |r) of the 
quantum system and is a LOMSO operator [24], and it is defined as D r = 
diag(0, 0, 1, 0, ...0), that is, the diagonal element (D r ) rr = 1 only for the 
index r and (D r ) tt = for any other index t ^ r. It has been proved [17] 
that when the selective phase-shift operation C s {9) of the marked state \s) 
acts on any quantum state Ix)^) a phase factor exp(— i9) is generated only 
when the quantum state \x) = \s), otherwise the quantum state (x)^) keeps 
unchanged, 

C,(e)\x)\S)=exp(-ie6 ax )\x)\S). (7) 
When the selective phase-shift operation C s {9) is applied to any superposition 

N-l 

= E a x\x)\S) of the quantum system, the unitary evolution process of 

x=0 

the quantum system is just the same as Eq.(5) of the oracle unitary operation 
U (9), 



6 



JV-1 

C s (9)\V)=exp(-t9)a s \s)\S) + £ a x \x)\S). (8) 

x=0,xj^s 

Therefore, the selective phase-shift operation C s (8) is really equivalent to the 
oracle unitary operation U (9) in the unitary transformation and one can use 
exactly the selective phase-shift operation C s {6) to replace the oracle unitary 
operation U to analyze the unitary evolution process of the quantum system 
under the action of the oracle unitary operation U (6). However, it must be 
pointed out that the selective phase-shift operation C s (6) is only applied to 
the work qubits / and independent of the auxiliary qubits S, while the oracle 
unitary operation U {6) is related to the whole quantum system. 

In order to describe conveniently the unitary evolution process of the 
quantum system under the selective phase-shift operations one may use the 
quantum state vector, i.e., the unity-number vector {a s k = ±1} [17] to repre- 
sent explicitly the selective phase-shift operation C s {6). In the quantum state 
vector {a s k } representation the diagonal operator D s of the marked state can 
be expressed as 

D s = (i^ + a\I lz ) (g)(i£ 2 + a s 2 I 2z ) <g) ... <g)(±£7 n + a s J nz ). (9) 
It has been shown that the quantum state vector {a s k } also determines uniquely 
the quantum state \s) [17]. Therefore, one may design a quantum search al- 
gorithm or sequence to determine directly the quantum state vector {a s k } 
and then obtains further the marked quantum state \s) instead of finding 
directly the marked quantum state \s) by a quantum search algorithm. This 
is different from the conventional quantum search strategy [3, 11]. This is 
also the starting point of the quantum search algorithms or sequences based 
on the quantum mechanical unitary dynamics. Obviously, the oracle unitary 
operation U o {0) contains the desired information to determine certainly the 
solution \s) to the quantum search problem. This also can be seen from 
the fact that the oracle unitary operation U (9) is equivalent to the selec- 
tive phase-shift operation C s (9), while the latter is determined uniquely by 
the quantum state vector {a s k }. Now an NMR ensemble quantum computer 
operating on the highly mixed states of a spin ensemble is used to solve the 
quantum search problem. As shown in Eqs.(2), (3), (5), and (8) the oracle 
unitary operations Uf and U {9) are independent of any quantum state \x) 
of the work qubits / but are related to the quantum state \S) of the auxiliary 
qubits S. When a quantum search algorithm based on the pure states and 
built up with the oracle unitary operations and oracle-independent unitary 
operations is used to determine the marked state \s) to the search problem 
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on a spin ensemble, the quantum state \S) of the auxiliary qubits S in the 
algorithm should be retained or replaced by the effective pure state in the 
spin ensemble [12], but any mixed states of the work qubits / of the spin 
ensemble can be used as the input states of the algorithm. For example, 
for the oracle unitary operation U o {0) the auxiliary quantum state |5')(S'| 
(] S) = |0)|1) (see Eq.(3)) in the matrix representation can be expressed as 

\S)(S\ = |0>(0|(g)|l>(l| = \E+\{S lz -S 2z )-S lz S 2z . (10) 
Then the auxiliary quantum state |5')(>S'| can be replaced by the effective pure 
state taking the same form as Eq.(10) and the initial quantum state \x) of 
the work qubits / of the algorithm can be replaced even by any highly mixed 
state of the work qubits / in the spin ensemble. However, the mixed state 
must be related to the marked state \s) in order that the quantum state vector 
{a s k } of the marked state |s) can be determined certainly by the algorithm. 
Note that any magnetization operator I kfl {k = 1,2, ...,n; p = x,y) and the 
diagonal operator D s of the marked state do not commute. It is easy to prove 
that the following density operators p(0) of the spin ensemble always contains 
the marked state \s), 

p(0) = Pl (0)\S)(S\ = (aE+e k I k ,)\S)(S\ (k = 1,2, ...,n; p = x,y), (11a) 
and 

n 

p(0) = Pl (0)\S)(S\ = (aE+ £ e k I k J\S)(S\ (116) 

fc=i 

where the operator E is unity operator, e k is the spin polarization parameter 
of the spin k in the spin ensemble, and a is the normalized constant. 

Now one can analyze conveniently the time evolution process of a quan- 
tum spin system or its spin ensemble under the action of the oracle unitary 
operation U o {0) by using the following general unitary transformation that 
describes the action of the selective phase-shift operation C s (8) on any den- 
sity operator P/(0) of the work qubits / of the spin ensemble [28]: 
CWpMCiO)- 1 = pM - (1 - cos0)[ P/ (O), D s } + 

+i sin 6[ Pl (0), D s ] + [(1 -cos #) 2 +sin 2 9}D sPl (0)D s (12) 
where the commutation [pj(0), D s ] + = pj(0)D s + D s pj(0) and [pj(0),D s ] = 
p j(0) D s — D s p j(0) . A general unitary transformation with a sequence of many 
selective phase-shift operations {C k (6 k )} can be derived directly from the 
basis unitary transformation (12), 

U o (0o, 6i, m -i)Pi(t)U o (6 , 9i, ^ m -i) _1 

m—1 m—1 

= Pi(t)-\p I (t), E (l-cos9 k )D k ] + + il Pl (t), E D k sm6 k ] 

k=0 k=0 
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m— ltn- 1 

+ E E K 1 " c os fc ) (1 - cos 6i) + smO k sm6 l ]D kPl {t) A 

fc=0 £=0 
m—1 

+ £ ^sin^(l-cos^)-sin^(l-cos^)]( J D fc p / (t)A-AP/(t) J Dfc) (13) 

l>k=0 

where the diagonal unitary operation U (9o, 9\, 9 m -i) is a sequence of the 
selective phase-shift operations: 

m—1 

U o (9 ,9i, 6 m -i) — n C k {9 k ). 

k=0 

On the basis of the unitary transformation (12) one can design a simple 
quantum search sequence to solve efficiently the quantum search problem on 
a spin ensemble. First, the spin ensemble is prepared at the initial mixed 
state of Eq.(llb). Then the oracle unitary operation U (9) is applied to 
the spin ensemble. The spin ensemble will evolve according to the unitary 
transformation (12). This unitary transformation is given explicitly below 
by inserting the diagonal operator D s of Eq.(9) in the quantum state vector 
{a s k } representation and the initial density operator of Eq.(llb) with the 
index /i — y into Eq.(12), 

n n n 

U (9) £ e k I ky U (9)- 1 =Z e k I ky - (l-cos0) E C^i + a{I lz ) <g> ... 

k=l k=l k=l 

®(!-Efc-i + al-ih-iz) ®{e k h y ) + a s k+1 I k+lz ) (g) ... 

n 

®{\E n + a s J nz ) -sm9 £ {\E X + a\I lz ) (g) ... (g)(|^_i + a s k _J k _ lz ) 

k=l 

®(e k a s k I kx ) ®{\E k+1 +a s k+1 I k+lz ) (g) ... ®{\E n +a s n I nz ) (14) 
where the term proportional to the unity operator E and the auxiliary quan- 
tum state \S) (S\ of the initial density operator of Eq.(llb) are omitted with- 
out losing generality. Note that here the auxiliary quantum state |S')(S'| is 
taken as the effective pure state and is of the LOMSO operator subspace, 
and hence it will keep unchanged when a z— direction gradient magnetic field 
is applied to the spin ensemble [6, 7]. Obviously, the generated mixed state 
of Eq.(14) under the action of the oracle unitary operation U a {9) on the spin 
ensemble is quite complicated, but it can be simplified by using the selective 
decoherence manipulation. First, a 90 degree hard pulse along y— direction 
is applied to all the spins of the work qubits /. Then a z— direction gradient 
magnetic field [6, 7] is applied to the spin ensemble to cancel all the nonzero- 
order multiple-quantum coherence including the single-quantum coherence, 
and then a zero-quantum dephasing pulse [29] is applied to the spin ensemble 
to cancel all the zero-quantum coherence.. In the final the density operator 
of Eq.(14) is reduced to the form 
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Uo{0) £ SkhyUoie)' 1 =>p f = (2/N) £ e k a s k I kz (15) 

k=l k=l 

Obviously, the density operator pj- of Eq.(15) carries the desired information, 
i.e., the quantum state vector {a|}, to solve the quantum search problem. 
Now, a read pulse 90° x is applied to all the spins of the work qubits / and 
then the NMR receiver is turned on to record in phase- sensitive mode the 
NMR signal of the density operator pj of Eq.(15). In order to cancel the effect 
of the auxiliary qubits S of Eq.(10) on the recorded NMR signal the spins 
of the auxiliary qubits S are decoupled during recording the NMR signal. 
From the fast Fourier transform spectrum of the recorded NMR signal one 
can determine certainly the quantum state vector {a s k } when the signal-to- 
noise ratio of the NMR signal of Eq.(15) is high enough. The density operator 
of Eq.(15) shows that the NMR signal intensity is proportional to the factor 
(2/N), indicating that the NMR signal intensity reduces exponentially as the 
qubit number of the work qubits I. This characteristic feature is similar to 
the classical NMR quantum computation based on the effective pure states 
[15]. Therefore, the above quantum search algorithm can only be used to 
solve efficiently the quantum search problem with a few qubits. 

However, the significant difference between the quantum search algo- 
rithms based on the spin ensemble and the pure quantum state [3] or the 
effective pure state [11] is that the Grover quantum search algorithm based 
on the pure quantum states [3] or the effective pure state [11] can achieve 
only quadratic speedup over the classical counterparts due to the limit of 
quantum mechanical measurement [6], which is a limit in quantum mechan- 
ical principle, that is, the conventional Grover algorithm is a quadratic but 
not polynomial-time algorithm, but the ensemble quantum search sequence 
can achieve an exponential-time speedup over the classical counterparts, that 
is, this sequence is a polynomial-time one (one performs only twice times the 
oracle unitary operation Uf in the algorithm) and really moves the frontier 
between the solvability and intractability of complexity problems. Though 
the final output NMR signal of the ensemble quantum search sequence re- 
duces exponentially as the size of the solved search problem, this limit is 
in techniques but not in quantum mechanical principles and therefore, this 
limit could be overcome as the relevant techniques are improved [7, 8, 22, 
23, 24]. There is also another advantage for the ensemble quantum search 
sequence that the sequence can be run even on a spin system with a short 
decoherence time since in the whole sequence the oracle unitary operation 
Uf is performed only twice times. 
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There may be two possible schemes [18] to improve the above ensemble 
quantum search sequence so that it could be used to solve the search problem 
with a larger number of qubits. One possible scheme is that the initial spin 
polarization of a spin ensemble is enhanced sufficiently with various possi- 
ble NMR signal enhancement techniques including spin polarization transfer 
techniques [6, 7, 22] and other techniques [23, 24] so that the NMR signal 
of the density operator of Eq.(15) still can be observable even in the spin 
system with a larger number of qubits, although the final NMR signal of the 
sequence is reduced exponentially as the qubit number of the spin system. 
This scheme retains the original advantage of the sequence that it still can 
be run on the spin systems with a short decoherence time. However, this 
scheme is rigorous in technique. Another scheme is that one could enhance 
the final detectable NMR signal by running a quantum search sequence that 
consists of a number of the oracle unitary operations and oracle-independent 
unitary operations on the spin ensemble. But this scheme requires that the 
spin ensemble must have a long decoherence time. Obviously, there is a max- 
imum call number M s of the oracle unitary operation U (9) or a maximum 
running time T 2max of the quantum search sequence in a spin ensemble due 
to the decoherence effect. The final output NMR signal of the sequence will 
decrease simply but not be enhanced again as increasing the call number of 
the oracle unitary operation in the sequence if the call number is larger than 
the maximum one M s . This is because in that case running the sequence 
will take more time than T 2max , so that the decoherence effect becomes the 
dominant effect negatively upon the final NMR signal. 

A simple method to enhance the output NMR signal by many calls of 
the oracle unitary operation is the conventional spin echo experiments [5, 
7]. The qubit interaction in the diagonal operator D s of the marked state 
\s) can be decoupled one by one by the spin-echo experiments. First, any 
Hamiltonian D s (n — k) (k = 1,2, ...,n — 1) that is independent of k qubits 
of the spin system is prepared with the spin echo sequences by starting the 
diagonal operator D s , 

D s (n - 1) = D s + exp(-inI nx )D s exp(inl nx ), 

D s (n - 2) = D s (n - 1) + exp(-i7r/ n _i x )D s (n - 1) exp(i7r/ n _i x ), 

D s (n-k) = D s (n-k + l) + exjp(-inI n - k+lx )D s (n-k + l)exp(i7rI n - k+ i x ), 

n—k 

where D s (n — k) =0 (\Ei + afli z ). Then one builds up a new auxiliary 

i=i 

oracle unitary operation V os {k,6) = exp[—i9D s (n — k)] with the Hamilto- 
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nian D s {n — k). There are 2 k oracle unitary operations U (9) in the unitary 
transformation V os (k, 9). Now using the oracle unitary operation V os (k, 9) to 
replace the oracle unitary operation U Q {9) in Eq.(14) to perform quantum 
search one can find that the final NMR signal intensity of Eq.(15) is propor- 
tional to the factor 2 fc , indicating that the output NMR signal intensity of 
the new quantum search sequence increases linearly as the call number of the 
oracle unitary operation. Obviously, this is an ineffective sequence to extend 
the realizable size of the NMR ensemble quantum computation. An improved 
approach to the simple spin-echo experiments is that following the action of 
the unitary operation V os {k,9) on the initial density operator one applies a 
sequence of non-oracle selective unitary operations Yl C r {9) to the spin en- 

semble. Then the NMR signal intensity of Eq.(15) is proportional to the sum 
S £ k{^2 °|+ S a i)Ikx, where the quantum state vectors {a s k } and {a r k } are 

k s r 

related to the oracle and the non-oracle selective unitary operations, respec- 
tively. Now, if the factors a s k and a k have the same sign then the final 

s r 

NMR signal is enhanced, otherwise it is reduced. Therefore, the NMR signal 
intensity may be enhanced when the non-oracle selective unitary operations 
are chosen properly. The more the number of the non-oracle selective uni- 
tary operations is in the sequence, the stronger the NMR signal is. However, 
there is a question that how many non-oracle selective unitary operations can 
be used to enhance effectively the NMR signal under the condition that the 
quantum state vector {a s k } can be extracted certainly from the final NMR 
signal intensity. This question limits really the maximum enhancement of 
the NMR signal by the sequence since the more the non-oracle selective uni- 
tary operations are used in the sequence it becomes more difficult for the 
solution information of the search problem to be extracted from the NMR 
signal. In order to improve further the method the NMR multiple- quantum 
spectroscopy is proposed to solve the search problem in the following. I will 
discuss the quantum search NMR multiple-quantum spectroscopy from the 
view point of the NMR experiments and leave its complexity problem in the 
future. 



3. The NMR multiple-quantum spectroscopy 

The basic idea for the multiple-quantum NMR spectroscopic method used 
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to perform the quantum search is that information of the unknown marked 
state \s) is first loaded on the multiple-quantum coherence of a spin en- 
semble by using iterative sequences of the oracle unitary operations U (6) 
and oracle-independent unitary operations, the search problem then could 
be solved by extracting the desired information of the marked state from 
these NMR multiple-quantum transition spectra obtained by NMR measure- 
ment techniques. The important point for the method is that the multiple- 
quantum coherence of the spin ensemble need first to be excited effectively 
by a sequence of the oracle unitary operations and oracle-independent uni- 
tary operations. Here a general (formal) quantum search sequence based on 
the multiple-quantum NMR spectroscopy [7, 8] is proposed. The sequence is 
constructed with the oracle unitary operations and oracle-independent uni- 
tary operations and may be divided into the following several steps: 

(a) The suitable initial input state for the quantum search sequence is pre- 
pared at the mixed state of a spin ensemble such as Eq.(lla) or (lib), 

n 

p(0)= Pl (0)\S)(S\ = (a E+ £ e k I kp )\S)(S\ (p = x,y,z). 

k=l 

(b) Apply the iterative sequence U(D S ) which is built up with the oracle 
unitary operations U {6) and oracle-independent unitary operations to the 
initial mixed state p(0) to excite the multiple-quantum coherence, 

p(n) = U(D s )p(0)U{D s ) + 

(c) During the t\ evolution period under the effective spin Hamiltonian H 
the multiple-quantum coherence are labelled by their own frequencies so that 
various order multiple-quantum coherence can be distinguished clearly in 
NMR multiple-quantum spectra, 

p(ti,Ti) = exp(-i#£i)p(Ti) exp(iifii) 

(d) Apply the second iterative sequence V(D S ) of the oracle unitary oper- 
ations and oracle-independent unitary operations to convert the multiple- 
quantum coherence into the single-quantum coherence so that they can be 
detected by NMR probes. 

p{T2,t 1 ,T 1 )=V(D 9 )p(t 1 ,T 1 )V(D a ) + 

(e) The NMR signal detection along with q— direction (q = x, y, z) 

n 

P f = P f (h) = Tr{J2 hqP^tiTi)}. 
k=i 

The longitudinal component l kz in the final NMR signal p^ is not observable 
in NMR experiments, but it becomes observable if an additional 90 degree 
pulse converts it into the single-quantum coherence. The multiple-quantum 
coherence are usually detected indirectly by NMR measurement techniques 
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such as the two-dimensional NMR method [7, 8]. The usual two-dimensional 
NMR experiment methods have been developed to do quantum computation 
on the basis of the concept of the effective pure states [12], but here NMR 
multiple-quantum spectroscopy is emphasized. The oracle iterative sequence 
U(D S ) carries the information of the marked state \s) and in the step (b) the 
information is really transferred into the multiple-quantum coherence of the 
spin ensemble through applying the iterative sequence U(D S ) to the initial 
density operator of the spin ensemble. The multiple-quantum coherence (or 
transitions) are detected indirectly by the two-dimensional NMR measure- 
ment techniques [7, 8], which is achieved by the steps (c), (d) and (e) in 
the sequence. Obviously, the output NMR signal Pf(t±) carries the informa- 
tion of the marked state \s). Then by fast Fourier transforming the output 
NMR signal Pf{t\) on the time variable t\ one obtains the multiple-quantum 
spectrum. Now the spectral parameters including the multiple-quantum res- 
onance frequencies, phase, and intensities of the multiple-quantum spectrum 
contain the information of the marked state \s). Then one could extract the 
quantum state vector {a s k } and hence find the marked state \s) from these 
spectral parameters by NMR data processing. 

3.1 The expression for the output NMR signal 

For convenient discussion, here consider a homonuclear spin ensemble 
where the spin polarization parameter for all the spins in the ensemble are 
the same (e = Then the initial density operator of work qubits / in the 

n 

step (a) is written briefly as p/(0) = eF p — e ^2 Ik p without losing generality 

k=l 

since the auxiliary quantum state |5')(S'|, the term proportional to the unity 
operator E, and the spin polarization factor keep unchanged during the 
evolution process of the quantum search sequence. From the steps (a)-(e) of 
the sequence the final NMR signal intensity of the sequence can be generally 
expressed as 

pjfa) = Tr{F q V{D s )eM-iHt 1 )U{D s )eF p U{D s ) + e-K V {iHt 1 )V{D s ) + } 
= eTr{V(D s ) + F g V(D s ) exp(-iHt 1 )U(D 8 )F p U(D 9 ) + exp^)} (16) 
where the invariance of the trace operation Tr to cyclic permutation has been 
introduced. In high-resolution NMR spectroscopy the spin Hamiltonian H 
of the spin system in the spin ensemble in a high magnetic field is a longitu- 
dinal magnetization and spin order (LOMSO) operator and usually can be 
written as [7, 8] 



14 



n n 

H = E Qkhz+ E ^Jkihzhz (17) 

k=l k>l 

With the help of the eigenbase \k) of the Hamiltonian H, H\k) — 0Jk\k) {h — 
1), which are also the usual quantum computational bases, the NMR signal 
of Eq.(16) can be expanded as 

Pf(h) = £ E Q*jk p jk exp(-iw ifc *i) (18) 

where the multiple-quantum transition frequency Uj k = Uj — uj k and the Her- 
mitian operators P(D S ) and Q(D S ) are defined as 

P{D S ) = U(D S )F P U(D S ) + , (19a) 
Q(D a )+ = Q{D S ) = V(D s )+F q V(D s ), (196) 
and the matrix elements Pj k and Qj k are given by 
P jk = {j\P(D 8 )\k),Q jk = {j\Q(D a )\k). 
Although the eigenvalue u k of the Hamiltonian H is independent of the quan- 
tum state vector {a|} of the marked state \s), the excited multiple-quantum 
transition frequencies {oojk} may contain the information of the marked state 
\s) besides the amplitudes and phases of the multiple-quantum transitions 
since the excitation of the multiple-quantum coherence by the oracle uni- 
tary sequence U(D S ) is dependent on the marked state. Actually, the oracle 
unitary sequence U(D S ) decides which orders multiple-quantum coherence 
are excited effectively and finally detected in the quantum search sequence. 
If the inphase multiple-quantum spectrum is required to be observed for a 
given quantum order one had better require that the oracle unitary sequence 
V(D S ) satisfies the relation [30]: 

Q(D s )+ = eM-i<pFz)P(D s )eMi<pFz)- (20) 
Then the NMR signal intensity of Eq.(18) can be further expressed as 

Pf(ti) = £ EE \ p jk\ 2 exp(irmp) exp(-tt^ fe ti) (21) 

m j,k 

where the quantum order m = Mj — M k and M k is the magnetic quantum 
number: F z \k) = M k \k) (H = 1). The NMR signal of Eq.(21) indicates that 
each m— order multiple-quantum transition has the same phase, although 
the multiple-quantum transitions with different quantum order m may have 
a phase difference. Therefore, by using separation techniques for different 
order multiple-quantum transitions such as TPPI technique [7, 8] one can 
record the inphase multiple-quantum peaks with their own resonance frequen- 
cies or quantum order in the multiple-quantum spectrum. The intensities, 
frequencies, and phases of the multiple-quantum peaks carry the informa- 
tion of the marked state. Then from these spectral parameters one could 



15 



extract the desired information of the marked state. Although the inphase 
multiple-quantum spectrum of Eq.(21) may be more convenient and useful 
in the measurement of structural properties of molecules in NMR multiple- 
quantum spectroscopy [30], the multiple-quantum spectrum of Eq.(18) may 
contain richer information of the marked state and may be more useful to 
solve the quantum search problem. 

If the effective spin Hamiltonian is prepared as H = uoF z instead of that 
one of Eq.(17) in the t\ evolution period, all the multiple-quantum coherence 
with the same quantum order have the same precession frequency. That 
is, for all the k— order multiple-quantum coherence their precession frequen- 
cies are the same and equal f k = kuj (k = 0, ±1, ±2, ±n). In this case 
there are (2n + 1) different order multiple-quantum resonance peaks with 
their own resonance frequencies {fk} in the multiple-quantum spectrum for 
an n— qubit spin system. Since different order multiple-quantum coherence 
in the density operator Pf(ti) of Eq.(18) are labelled with their own tran- 
sition frequencies (or precession frequencies), they can be distinguished in 
the NMR multiple-quantum frequency spectrum. It is known that num- 
ber of independent multiple-quantum coherence increases exponentially as 
the qubit number in a spin system [7]. If each multiple-quantum coherence 
has its own precession frequency in a spin system, then amplitude of each 
multiple-quantum coherence will decrease exponentially as the qubit num- 
ber and this will make each multiple-quantum spectral peak unobservable 
in the multiple-quantum spectrum of a spin system with many qubits since 
the total amplitude of the multiple-quantum coherence of a spin system is 
fixed and is determined by the initial density operator. The importance to 
choose the spin Hamiltonian H = ojF z to label the multiple-quantum coher- 
ence in the t\ evolution period instead of the spin Hamiltonian of Eq.(17) 
therefore becomes very clear: all the same order multiple-quantum coherence 
have the same precession frequency and their resonance frequency peaks in 
the multiple-quantum spectrum can be coherently co-added in amplitude and 
become a single spectral peak, as can be seen in the NMR signal of Eq.(18), so 
that the multiple-quantum spectral peak may increase greatly in a favorable 
case. This shows that at least some lower order multiple-quantum spectral 
peaks in the spin ensemble do not reduce exponentially as the qubit number 
if the initial density operator is excited effectively to the multiple-quantum 
coherence by the oracle iterative sequence U(D S ). Therefore, complexity of 
the quantum search sequence now is only dependent on the excitation of the 
multiple-quantum coherence by the oracle iterative sequence U(D S ) and the 
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extraction of the solution information of the search problem by the NMR 
data processing. Of course, it could become more difficult to extract the 
information of solution of the search problem using this method labelling 
multiple-quantum coherence with respect to that one using the Hamiltonian 
of Eq.(17). By the NMR signal expression (18) one can explicitly calcu- 
late the output NMR signal intensities of various order multiple-quantum 
transitions, 

Pf(.h) = E [e E' Q* jk Pjk\ expi-imuh) (22) 

m=—n j,k 

where the sum E runs on ly those indexes (J, k) with the quantum order 

m = Mj — M k . The m— order multiple-quantum peak intensity and phase 

therefore are determined from 

I{m) = [e E' Q* jk Pjk}. (23) 
j,k 

Now the information of solution of the search problem could be extracted 
from the spectral parameters I{m) by NMR data processing. 

3.2 The multiple-quantum operator algebra spaces and the con- 
struction of the oracle unitary sequences 

The multiple-quantum operator algebra space [24, 25] states that the 
whole Liouville operator space of a spin ensemble contains the even-order 
multiple-quantum operator subspace, the zero-quantum operator subspace, 
and the longitudinal magnetization and spin order (LOMSO) operator sub- 
space and each former subspace contains the followed subspaces. For ex- 
ample, the zero-quantum operator subspace contains the whole LOMSO 
operator subspace. There are important properties for the zero-quantum op- 
erators and the even-order multiple-quantum operators. Any operator of the 
Liouville operator space will keep its quantum order unchanged when it is 
acted on by a zero-quantum operator and the parity of quantum order of an 
operator of the Liouville operator space can not be changed when the oper- 
ator is acted on by an even-order multiple-quantum operator [24, 31]. With 
the help of the closed properties of these operator subspaces one may sim- 
plify the construction of the oracle unitary sequences U{D S ) and V(D S ). In 
general, the oracle unitary sequences U(D S ) and V(D S ) are composed of the 
oracle unitary operation U {6) and oracle-independent unitary operations. 
Suppose that the oracle unitary sequence U(D S ) is constructed in the form: 
U(D S ) = R k O k - 1 R k ^ 1 ...0 1 R 1 (24) 
where each unitary operation R k is an oracle-independent unitary operation 
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and 0^ = Ok{D s ) an oracle unitary operation sequence which contains the 
oracle unitary operation U Q (9). Now assume that the initial mixed state Pj{0) 
in the quantum search sequence is taken as the longitudinal magnetization 

n 

operator, p 7 (0) =Y^ ^khz, which is also a zero-quantum operator. There 
k=i 

may be two simpler schemes to load the information of the marked state 
\s) on the multiple-quantum spectrum. One scheme is that only the zero- 
quantum transitions are used to carry the information of the marked state. 
The zero-quantum operator subspace is the smaller operator subspace in 
the Liouville operator space. Then it is a better scheme to load the infor- 
mation of the marked state on the zero-quantum operator subspace since 
this could simplify the loading and extraction of the desired information of 
the marked state. Obviously, this requires that the oracle unitary sequence 
U(D S ) should be prepared as a zero-quantum unitary operator. When the 
oracle unitary sequence U(D S ) is applied to the initial mixed state p 7 (0) 
only the zero-quantum coherence of the spin system can be excited accord- 
ing to closed property of the multiple-quantum operator algebra space since 
the initial mixed state P/(0) is also a zero-quantum operator. In order to 
distinguish the effects between the oracle unitary operations {Ok} and the 
oracle-independent unitary operations {Rk} on the excitation of the zero- 
quantum coherence it could be better to choose the unitary operations {Rk} 
as the LOMSO unitary operators and the oracle unitary operations {Ok} 
as the zero-quantum unitary operators. It is known from multiple-quantum 
operator algebra space formalism that any pair of the LOMSO operators 
commute with each other, and the closed property of the zero-quantum op- 
erator subspace shows that any zero quantum operator is transformed uni- 
tarily into a zero-quantum operator when it is acted on by any zero-quantum 
unitary operator. Then the zero-quantum coherence will be excited only by 
the oracle unitary operations {Ok}, while the LOMSO unitary operations 
{Rk} should have not net contribution to the excitation of the zero-quantum 
coherence.. But the unitary operators {Rk} could be chosen properly so that 
the excitation efficiency of the zero-quantum coherence can be enhanced. If 
the zero-quantum coherence spectra are observed in the output NMR sig- 
nal of the quantum search sequence it can be sure that the zero-quantum 
spectra contain the information of the marked state. Another scheme is that 
the nonzero order multiple-quantum coherence and especially the even-order 
multiple-quantum coherence are used to carry the information of the marked 
state. According to the multiple-quantum operator algebra space formalism 
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all the even-order multiple-quantum operators form a closed operator sub- 
space of the Liouville operator space. The closed property of the even order 
multiple-quantum operator subspace shows that the initial density operator 
(a LOMSO operator) will be transferred into even order multiple-quantum 
coherence if it is acted on by an even-order multiple-quantum unitary oper- 
ator. Therefore, if all the oracle unitary operations {Ok} are taken as the 
even-order multiple-quantum unitary operators and the oracle-independent 
unitary operations {Rk} as the proper LOMSO or zero-quantum unitary 
operators, then even order multiple-quantum coherence will be excited by 
the oracle unitary sequence U(D S ) in the spin ensemble and moreover, the 
nonzero even-order multiple-quantum coherence can only be excited by the 
oracle unitary operations {Ok} and the oracle-independent unitary opera- 
tions {Rk} will have not net contribution to the excitation of the nonzero- 
order multiple-quantum coherence when the oracle unitary sequence U(D S ) 
is applied to the initial density operator p/(0). Therefore, the information of 
the marked state now is contained mainly in the nonzero even order multiple- 
quantum spectra of the output NMR signal of the sequence. How easy to ex- 
tract the information {a s k } of the marked state from these multiple-quantum 
spectra? This depends on what the explicit forms of the oracle-independent 
operators {Rk} and the oracle unitary operations {Ok} are taken. The im- 
portance is that the oracle unitary sequence U (D s ) needs to be constructed 
properly so that the multiple-quantum coherence can be excited effectively. 

(a) the preparation for the zero-quantum unitary operators 

There are a number of methods to construct the zero-quantum unitary 
operations with a sequence of the oracle unitary operations U {6) and oracle- 
independent unitary operations. One simple and systematic method is that 
the preparation can be performed by the discrete Fourier analysis and the 
phase cycling technique [17] by starting the diagonal operator D s of the 
marked state \s). According to the multiple-quantum operator algebra space 
any diagonal operator Di, including the diagonal operator D s of the marked 
state \s), can be expressed as a linear combination of the LOMSO basis 
product operators {Z k } = {E, I kz , 2I kz I lz , 4:I kz I lz I mz , 2 n - 1 I lz I 2z ...I n z} 
including the unity operator E (E — Z ) [24], 

N-l 

D k = E (A) kl Z t . (25) 
i=o 

Then the Hermitian product operators X\ = exp(— i^F y )Zi exp(i^F y ) can be 
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constructed by 

N-l 

X l = 2 l ~ l I klX I k2X ...I klX = exp(-i§ F y )[ £ (A-^jDj] exp(if F„) (26) 

and in a more general case any linear operator function f(D s ) related to 
the diagonal operator D s of the marked state \s) and any other diagonal 
operators D r by 

/(£>.) =EW 
z=i 

JV-1 

= exp(-if b t £ (A" 1 )^] exp(if F,) (27) 

1=1 j=0 

where b\ is a coefficient. Now by starting the operator function of Eq.(27) one 
can prepare the Hermitian zero-quantum operators with aid of the discrete 
Fourier analysis and the phase cycling technique [17], 

iVi-l 

H zq (D s ) = ±- £ exp(-iip k F z )f(D 9 )exp(-i<p k F g ) 

k=0 

l>k=l 
q>p>l>k = l 

+J« M (K})(^r W + vtW) 

+^ pq ({a%})(l£lrW+ItfW} + ---> (28) 
where the phase <y2 fc = k2ir/Ni, k = 0, 1, 2, iVi — 1 and the integer iVi > n. 

By the N\ — step phase cycling systematically only the zero-quantum opera- 
tor H zq (D s ) of the operator function f(D s ) is retained but all nonzero order 
multiple-quantum operators are cancelled. The zero-quantum unitary op- 
erator then can be constructed generally with the Hermitian zero-quantum 
operator of Eq.(28), 

U zq (D s , t) = exp{-iH zq (D s )t]. (29) 
Since the operator function f(D s ) contains the diagonal operator D s of the 
marked state \s) the zero-quantum operator H zq (D s ) of Eq.(28) carries the 
information of the marked state \s). The zero-quantum unitary operator of 
Eq.(29) may be further decomposed into a sequence of the oracle unitary 
operations and oracle-independent unitary operations by using the Trotter- 
Suzuki formula [32, 33]: 

lim ( e - iA i/ m e - iA 2/ m e -iA n /m^m _ e ~i(A 1 +A 2 +...+A n ) (30) 
m— »oo 

Now the decomposition (30) of the zero-quantum unitary operator of Eq.(29) 
may be used as the oracle unitary sequence U (D s ) to excite the zero-quantum 
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coherence in the quantum search sequence. 



(b) the preparation for the nonzero-order multiple-quantum uni- 
tary operators 

Using the discrete Fourier analysis and the phase cycling technique [17] 
by starting the diagonal operator D s of the marked state \s) one still can 
construct generally any nonzero order multiple-quantum coherence operators 
and further uses the operators to construct the multiple-quantum unitary 
operators. But a simple scheme to prepare a nonzero-order multiple-quantum 
unitary operator as a sequence of the oracle unitary operations and oracle- 
independent unitary operations is given below. For any Hermitian product 
operator Xi = 2 l ~ 1 I klx I k2X ...I klX of the x— direction product operator set {E, 
I kx , ^Ikxhx-, ^Ikxhxlmx, • ••} one may construct a number of the nonzero order 
multiple-quantum operators such as i[Xi, (D ±Djv-i)]±- The typical two of 
which are given explicitly by 

i[X h (D - D N _ t )] 
i 

= +\i (/*,* - iI kjV ) ®{\E kl+1 + h l+1 z) - ®{\E kn + I knZ ) 
3=1 

i 

-\% (h jX + ih iV ) ®{\E kl+1 - h l+lZ ) ... ®{\E kn - I knZ ) 

3=1 
I 

-\i (h jX + ih jV ) ®{\E kl+1 + h l+1 z) - 0(|£fc„ + h nZ ) 

3=1 

+\i {h jX -ih jy )®{\E kl+1 -I kl+lZ )®...®{\E kn -I knZ ), (31) 

3=1 

[X h (D + £iv-i)] + 
i 

= +\ (h jX - U kj y) 0(§^ +1 + h l+1 z) - <8>(lE kn + I knZ ) 
3=1 

+ \ (h jX + U kj y) ®{\E kl+1 ~ h l+1 z) - ®{\E kn - h n z) 
3=1 

+\ (h jX + U kj y) ®{\E kl+1 + h l+lZ ) - ®(\E kn + I knZ ) 

3=1 

+\ (h ]X - U kj y) ®{\E kl+1 - h l+lZ ) ... ®{\E kn - I knZ ). (32) 

3=1 

Obviously, the operators [X;, (D — -Djv-i)]± are /—order quantum opera- 
tors. Therefore, if the integer index I is an even number the operators 
[X ; , (Dq — D N -i)]± will be even-order multiple-quantum operators. Note 



21 



that when the integer index / is an even number the Hermitian operator 
Xi = 2 l ~ 1 Ik lx Ik 2 x---Ik l x also is really an even-order multiple-quantum oper- 
ator which is a linear combination of the even-order multiple-quantum op- 
erators with various quantum orders smaller than I including zero-quantum 
operators. In a more general case an /—order multiple-quantum operator Hi 
could be further generated from the Hermitian operators i[X h (D ± D^-i)] 
and [Xi, (D ± D N _i)] + , for example, 

Hi =E <*t[iXi, (D ± D N -!)}+ E PftXi, (D ± Av-i)]+, 

i i 
Hi = exp(-mI kl+]X )[iXi, (D ± Av-i)] exp(inI kl+jX ) (j = 1,2, ...,n- I). 

To load the information of the marked state \s) on the multiple-quantum 
operators one may first prepare the operators {X{\ or the operator func- 
tion f(D s ) from the diagonal operator D s of the marked state \s) and other 
diagonal operators by Eq.(26) or (27). Then these operators are further 
used to construct the desired multiple-quantum operators in a similar way 
to [Xi, (Dq — Djv-i)]±- Once the desired Hermitian multiple-quantum oper- 
ator Hi is prepared one can further construct the multiple-quantum unitary 
operator Ui(D s , A;) = exp(— iXiHi). There is a well-known operator identical 
relation for any two Hermitian operators A and B [33, 34]: 

lim (^ e iA /\/™ e iB /\/™ e - iA /Vm e -'i B /\/™^ m — g-[A-B] (33) 

This formula and Eq.(30) could be helpful for the decomposition of the 
multiple-quantum unitary operator into a sequence of the oracle unitary oper- 
ations and oracle-independent unitary operations. This decomposition could 
be used as the oracle iterative sequence U(D S ) in the quantum search se- 
quence. 

3.3 The Grover iterative sequence in the spin ensemble 

According the above scheme of Eq.(24) to construct the oracle unitary 
sequences U(D S ) a simple oracle unitary sequence U(D S ) is constructed ex- 
plicitly below. This sequence is similar to the Grover iterative sequence [3, 
17]. The time evolution propagator for the oracle unitary sequence (here 
denote U(D S ) as U(m)) is constructed simply by 

U(m) = [exp(-iiD I4 ) exp(-iirD x s )] m (34) 
where the oracle-independent unitary operation exp(— inD N _i) is a LOMSO 
unitary operator and the oracle unitary operation sequence exp(— inD*) a 
multiple-quantum unitary operator. The Hermitian multiple-quantum oper- 
ator D x a is defined as 
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D x s = exp{-t^F y )D s exp(t^F y ) = <g) {\E k +a%I kx ). (35) 

k=l 

By using the relation between the two diagonal operators D s and D [17], 
D s = exp(^F x )U ox (-^)D U ox ^)exp(t^F x ) (36) 

n 

where the auxiliary oracle unitary operation U ox {9) = Y[ exp[— i6a s k I kx ], one 

k=i 

can re-express the propagator of Eq. (34) as 
U(m) = exp(-qF y )exp(-qF x )U ox (-%) 

x [exp(-iirD$) exp(-inD )] m U ox (i ) exp(i§ F x ) exp(i^F y ). (37) 
It follows from the definition of the diagonal operator D s of the marked state 
that there is an expansion for the unitary operator exp(±inD s ) [17], 

exp(±i9D s ) = E+(-l+exp(±i9))D s . (38) 
The iterative part of the propagator U(m) of Eq.(37) then can be simplified 
by the expansion 

G(m) = [exp(-iTrDg) exp(-inD )] m 

= [E- 2D - 2D% + 4£>jJ£> ] m (39) 
It can turn out based on the definition of the diagonal operator D k that 
the operator set {E, D , D%, D Dq, DqDq} form a closed operator algebra 
subspace. Then the iterative sequence of Eq.(39) can be generally written as 
G{m) =E + ai(m)D + u 2 {m)D* + a 3 (m)D D% + a 4 (m)D%D . (40) 
Using the relations G(m + 1) = G(l)G(m) = G(m)G(l) one can obtain the 
recursion equations for the coefficients {a k (m)}, 

ai(m + 1) \ _ / -1 + A/N 2/N \ ( a x {m) \ ( -2 



a 3 (m+l) J V - 2 - 1 J V «3(m) 7 ' V 7' 

a 2 (m+ 1) \ _ / -1 -2/7V ^ f a 2 (m) ^ + / -2 \ (41fe) 



a 4 (m+l) y \2 -1+4/N J \ a 4 (m 
and cti(m) = a2(m) and azt(rn) = — {2a\{m) + a 3 (m)). The recursion equa- 
tions (41a) and (41b) have the following solution: 

N 

ai(m) = a 2 ( y m) = — — — -[1 — cos(m#)], 

a 3 (m) = — ^ -[— 1 + cos(m#) + ^/N — lsin(m#)], 

0:4(777.) = — -[1 — cos(m#) + \f N — 1 sin(m#)] 

where cos6» = -1 + 2/JV and sin6» = 2ViV - 1/N. 

The time evolution of a spin ensemble under the propagator of Eq.(37) 
then can be expressed as 

Pl (m) = U{m)p I {Qi)U{m)- 1 
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= exp(-if F y ) exp(-if F x )U ox (-% )G{m) 

x */«,(§) ex V {i\F x ) exp(qF y ) Pl (0) exp{-t^F y ) exp(-i§ F x )U ox {~l ) 
xG(m)+U ox (\) exp(*f F x ) exp(i|F J/ ). (42) 

n 

Now assume that the initial density operator is given by p/(0) =52 £khz- 

k 

Then one has the unitary transformation relations: 

U„{\) exp(if F,) exp(if F„)p,(0) exp(-i|F„) exp(-i§ jyir„(-§) 

fc 

n n 

G(m){- 52 e k I kx }G{m) + = - 52 £ k{hx + 7i(™)A) + 7 2 ( m )^o 

k k 

+<y 3 (m)D D% + <y 4 (m)D%D + >y 5 (m)D I kx + <y 6 (m)I kx D 
+l 7 (m)D%D I kx + -f s (m)I kx D D*} 
where the coefficients 7;(m) (/ = f , 2, 8) are given by 

7 1 (m) = 2^[ai(m)*a 3 (m) + Oi(to)o 3 (to)* + a 3 (m)*a 3 (m)}, 
7 2 (m) = |[o 2 (m) + o 2 (to)* + a 2 (^)*«2(^) 
+-^o 2 (m)*04(m) + ^o 2 (m)o 4 (m)*], 

7 3 (m) = \[a 3 (m) + oi(m)o 2 (m)* + o 2 (m)*a: 3 (m) + ^a 3 (m) 0:4(771)*], 
7 4 (m) = |[o 3 (m)* + oi(m)*o 2 (m) + o 2 (m)o 3 (m)* + ^a 3 (m)*04(m)], 
7s ( m ) = ai (m),7 6 (m) = oi(m)*, 
7 7 (m) = o 4 (m),7 8 (m) = o 4 (m)*. 
The density operator pj(m) of Eq.(42) then can be further given by 

n 

Pi( m ) = 52 £k{hz - li{m)D x s - 7 2 (m)D Ar _ 1 

k=l 

- 73 (m)L>f^_ 1 - 74 (m)£) JV _ 1 ^ + i b {m)D*I kz 

+ l6 (m)I kz D* + ll (m)D N _ l D x s I kz + 78 (m)/ te £)f (43) 

n 

Now the conversion efficiency of the initial density operator p 7 (0) =52 £ khz 

k 

to the multiple-quantum coherence under the Grover iterative sequence of 
Eq.(34) can be obtained from Eq.(43). The LOMSO operator in the density 
operator Pj{m) of Eq.(43) which have not any contribution to the multiple- 
quantum transitions is given by 

n 

Pz( m ) = 52 £ k{Ikz - jjl 1 (m)E - 7 2 (m)£)jv_i - ^ 3 (m)D N _ 1 
k=l 

-|74( m )^w-i+|75(m)^+^76( m )^+2]v77(m)^-i+^78( m )' D A'-i}' 
The density operator pj(m) is traceless since the initial density operator pj(0) 

is traceless. Then the LOMSO operator p z (m) is reduced to the form 



24 



p z (m) =E e k {I kz [l+^{ 1 ^m)+ 1& (m))]-^ ll (m)E+ ll (m)D N _ 1 }. (44) 

k=l 

In general, the diagonal operator D s defined by Eq.(9) can be expanded as 
a sum of the LOMSO product operator base 

n n 

D s = jf{E+ E (a 8 k 2I kz )+ E (a s k 2I kz )(a s 2I lz ) 

k=l l>k=l 

+ E (a s k 2I kz )(af2I lz )(a s m 2I mz ) 

m>l>k=l 

+ + (al2I lz )(a s 2 2I 2z ) «2/ n2 )}. (45) 

Then inserting the expansion (45) of the diagonal operator D N _i into Eq.(44) 
one obtains 

n 

Pz( m ) =E £klkz[l + jfilbi™) +7eM)] 
k=i 

Hit e*)7i("»){E (-24,)+ E (-2I kz )(-2I lz ) 

k=l k=l l>k=l 

n 

+ E (-2I kz )(-2I lz )(-2I mz ) 

m>l>k=l 

+ + (-2I lz )(-2I 2z ) (-2I nz )}. (46) 

Therefore, the conversion coefficient C m for the initial density operator (or 

n 

the initial longitudinal magnetization) pj(0) =E £ khz is given by 

k 

n 

C m = l + ±h 5 (m) + 7 6 M] - f(E £i/e fc )7iH- (47) 

The conversion coefficient C m = 1 means that the initial magnetization p/(0) 
is not transferred into any other components including multiple-quantum 
coherence in the spin ensemble, while C m = means that the initial mag- 
netization p/(0) is completely converted into other components in the spin 
ensemble. It can be seen from the coefficients {a k {m)} that the maximum 
7 1 (m), 7 5 (m), 7 6 (m) ~ 0(1) and occur at m ~ v^V- Then the conversion 
coefficient C m of Eq.(47) shows that the initial magnetization pj(0) can not 
be transferred effectively into multiple-quantum coherence for a large qubit 
number (N = 2 n ) and the transferred efficiency (I — C m ) decreases exponen- 
tially as the qubit number for any iterative number m. This result shows 
that the Grover iterative sequence (34) is not an effective sequence to en- 
hance the NMR signal and amplify the realizable size of a spin ensemble in 
solving quantum search problem in a spin ensemble. This is quite different 
from the pure-state quantum search algorithm which can achieve a quadratic 
amplification for the amplitude of the marked quantum state. One sees again 



25 



the difference between the ensemble quantum search computation and the 
pure-state quantum search computation. 

3.4 NMR signal enhancement by cross interaction between oracle- 
dependent and oracle-independent interactions 

The extraction of the solution information of the search problem in the 
multiple-quantum spectrum could be simplified by preparing the suitable or- 
acle iterative sequence U(D S ) and its effective Hamiltonian H(D S ) used to 
excite the multiple-quantum coherence.. As an example, assume that the 
information of the marked state \s) is loaded on the zero-quantum operator 
subspace of the Liouville operator space of the spin system. The preparation 
for the zero-quantum Hamiltonian carrying the information of the marked 
state can be carried out by using the discrete Fourier analysis and the phase 
cycling technique [17], as can be seen in Eq.(28). By starting the proper oper- 
ator function f(D s ) of the diagonal operator D s of the marked state and the 
oracle-independent operator function f(D r ) the zero-quantum Hamiltonian 
may be constructed generally 

iVi-l 

H ^ = W E eM-mF z )[f(D s ) + f(D r )]eM-mF z ) 

= H zq (D s ) + H zq (D r ) (48) 
and the corresponding zero-quantum unitary operator is already given by 
Eq.(29). If the oracle-independent Hamiltonian H zq (D r ) is suitably cho- 
sen, for example, it may be chosen as those zero-quantum coherence of the 
k— qubit (k < n) subsystem of the spin system and moreover, as the dominat- 
ing term in Eq.(48) with respect to the Hamiltonian H zq (D s ), then extraction 
of the solution information could be simplified and the zero-quantum coher- 
ence which carry the information could be enhanced in amplitude, as can 
be seen below. Now there are three types of the zero-quantum coherence 
in the spin ensemble. One type is of the k— qubit subsystem of the spin 
system, another of the rest (n — A;)— qubit subsystem, and the third is those 
cross zero-quantum coherence between the two subsystem. The last type 
is more important. Assume that all spins of the k— qubit subsystem have 
the same precession frequency and so do all spins of the rest (n — k)— qubit 
subsystem, but the two precession frequencies are different each other. Then 
the cross zero-quantum coherence between the two subsystems may have 
precession frequencies different from zero frequency that is the precession 
frequency of the zero-quantum coherence of each subsystem. Therefore, the 
cross zero-quantum coherence with nonzero precession frequencies and those 
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zero-quantum coherence of the two subsystems can be discriminated in the 
zero-quantum spectrum. Since these cross zero-quantum peaks also contain 
the information of the marked state one could obtain the information of the 
marked state from the spectral parameters of the cross zero-quantum peaks. 
The zero-quantum peaks of the k— qubit subsystem are strongest since the 
oracle- independent Hamiltonian H zq (D r ) is the dominating term in the to- 
tal Hamiltonian H zq , while those peaks of the (n — A;)— qubit subsystem are 
usually weakest and even can not be observable if a small call number of the 
oracle unitary operation is used in the sequence, but the information of the 
marked state may not be obtained from these zero-quantum peaks because 
of overlapping of the strong zero-quantum peaks without the information. 
However, amplitude for the cross zero-quantum coherence is dependent on 
both the oracle- dependent and oracle-independent Hamiltonians. Then the 
cross zero-quantum peaks could be much stronger than those zero-quantum 
peaks of the (n — k)— qubit subsystem and could be observable even when the 
peaks of the (n — A;)— qubit subsystem are too weak to be observable. There- 
fore, one may use these cross zero-quantum peaks with nonzero precession 
frequencies to determine the solution of the search problem. Obviously, this 
NMR signal enhancement is achieved by the cross interaction between the 
oracle-dependent and oracle-independent Hamiltonians. To see more clearly 
the enhancement mechanism one first expresses the zero-quantum unitary 
transformation in the interaction frame as 

U zq (t) = exp(-iH gq t) = exp[-itH zq (D r )]TeM-i I H' zq (D s ,t')dt') (49) 

o 

where T is Dyson time-ordering operator and the oracle-dependent Hamilto- 
nian in the interaction frame is given by 

H' zq (D s ,t) = exp[ttH zq (D r )]H zq (D s ) exp[-itH zq (D r )]. (50) 
Equation (49) may be used to build up the quantum circuit U zq (t) which 
equals U(D S ) in the quantum search sequence [37]. The interaction Hamilto- 
nian (50) is responsible for generating the cross zero-quantum peaks. It can 
be expanded formally as 

H' zq (D s ,t) = H zq (D s )+tt[H zq (D r ), H zq (D s )] 

-\t 2 [H zq {D r ), [H zq (D r ), H zq {D s )\] + .... (51) 
One sees that there are the cross interaction terms between the two Hamil- 
tonians H zq (D r ) and H zq (D s ) in the higher-order terms of the expansion 
except the term H zq (D s ). These cross interaction terms are responsible for 
generating the cross zero-quantum peaks. Obviously, the cross interaction in 
strength is dependent on the oracle independent Hamiltonian H zq (D r ), and a 
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strong Hamiltonian H zq (D r ) will generate a strong cross interaction. There- 
fore, a strong oracle-independent Hamiltonian would enhance greatly the 
cross zero-quantum peaks since strong cross interaction would create strong 
cross zero-quantum peaks. On the other hand, the oracle-independent Hamil- 
tonian H zq (D r ) becomes more and more important in the higher order terms, 
as can be seen in Eq.(51). A cross interaction term with a high order such 
as [H zq (D r ), [H zq (D r ), H zq (D s )]...] could amplify the effect of the oracle- 
independent Hamiltonian on the cross zero-quantum peaks. 

Generally, the Baker- Campbell-Hausdorff (BCH) formula [38] can be ex- 
ploited to construct the higher-order cross interaction. For convenience, op- 
erators A = H(D r ) and B = H(D S ) denote the oracle-independent and 
oracle-dependent Hamiltonians, respectively. The operator A may be any 
Hermitian operator of the spin system and usually is taken conveniently as 
the LOMSO operator. The oracle-dependent Hamiltonian H(D S ) is usually 
taken as an Hermitian multiple-quantum operator. The BCH formula [35, 
36, 38] shows that the symmetric composition for the two operators A and 
B with the second order approximation is written as 

Sa{x) = exp(xA') exp(xB) exp(xA') = exp(xAx) (52a) 
with the Hermitian operator, 

xA l = x(A + B) + lx 3 ([B, [B, A']] - [A', [A', B]]) + (526) 
where x — it is the imaging time and A' = A/2 (the same denoting for 
operator B', A[, below). There is also another second order symmetric 
composition with a different product order of the operators A and B, 

Sb(x) = exp(xB') exp(xA) exp(xB') = exp(xBi) (53a) 
with the Hermitian operator, 

xB 1 = x(A + B) + \x\[A, [A, B']] - [B 1 , [B', A]]) + .... (536) 
Because the unitary operators S q (x) (q = A, B) of Eqs.(52a) and (53a) satisfy 
the time symmetry: S q (x)S q (— x) — E (E is unity operator) there are not 
even-order commutators such as [A, B], [A, [A, [A, B}}}, etc., in the composed 
operators A± and B x of Eqs.(52b) and (53b) [36]. Here an n— order cross 
interaction term is defined as a (n + 1)— order commutator of the operators 
A and B. Then the second-order cross interaction, which is defined as the 
third order commutators such as [A, [A,B]] and [B, [B, A]], can be built up 
in the same symmetric form as Eq.(52a) with the operators Ai and B\ again 
according to the BCH formula, 

exp(x^4 / 1 ) exp(-xBi) exp(xA' 1 ) = exp(xA 3 ) (54a) 
where the cross interaction A 3 is given by 

xA 3 = x(A 1 - B,) + lx s ([-B u [-B u A[]] - [A[, [A' v -B,}}) + ... 
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= lx 3 ([B,[B,A]]-[A,[A,B]]) + ..., (546) 

and 

exp(xB[) exp(— xAi) exp(xB[) = exp(xB 3 ) (55a) 
where the cross interaction B 3 is given by 

xB 3 = x(-A 1 + B x ) + |x 3 ([-A, [-A x , B[]] - [B[, [B[, -A,]]) + ... 

= lx*([A,[A,B]}-[B,[B,A]]) + .... ^ (556) 

It can be seen in Eqs. (54b) and (55b) that the operators A and B disappear 
and the dominating terms become the third order commutators [A, [A, B}} 
and [B, [B,A\] in the cross interaction A 3 and B 3 . If those commutators 
with order higher than four can be neglected in Eqs. (54b) and (55b) then the 
operators A 3 and B 3 are approximated as the second order cross interaction. 
The fourth-order cross interaction can be prepared by starting the second 
order cross interaction A 3 and B 3 , 

exp(xA' 3 ) exp(xB 3 ) exp(xA' 3 ) = exp(xA 5 ), (56a) 

and 

exp(xB' 3 ) exp(xA 3 ) exp(xB' 3 ) = exp(xB 5 ). (566) 
The higher order cross interaction also can be easily built up with the same 
method above. The complexity of the quantum circuit to construct the higher 
order cross interaction is usually dependent on the two Hamiltonians H(D r ) 
and H(D S ). However, it is easy to see that preparing 2m— order cross in- 
teraction terms A 2m +i and B 2m +i need |(3 m+1 — 1) and |(3 m+1 + 1) oracle 
unitary sequences exp(x.B), respectively. How high order cross interaction 
^2m+i (B 2m +i) is required for an n— qubit spin system so that the initial den- 

n 

sity operator p/(0) =^2 e^l^ z can be transferred efficiently into the multiple- 

k 

quantum coherence by the oracle unitary sequence U(D S ) = exp(— itA2 m +i) 
(exp(— iti?2m+i))? This is also dependent on the two Hamiltonians H(D r ) 
and H(D S ) but it could be lower than n for the proper H(D r ) and H(D S ). A 
much more efficient method to prepare the higher order cross interaction can 
be seen in [35, 36]. According to the method [35, 36] one first composes in 
a symmetric form the (2m — 1)— order composed unitary operators f2m-A x ) 
and fim-i^) by using the second order composed unitary operators Sa(x) 
and Sb{x) of Eq.(52a) and (53a), respectively, 

fL-M = exp[x(A + B)+ x 2m+1 A 2m+1 + ...] 

= S A (p 1 x)S A (p2x)...S A {Prx), (57a) 

fL-i{x) = exp[x(A + B) + x 2m+1 B 2rn+1 + ...] 

= S B (p 1 x)S B (p2x)...S B (Prx) (576) 



29 



where the parameters {p r } must be determined suitably and satisfy Pk = 

1, pj = Pr+i-j (j = 1, 2, r — 1), and it is usually not easy to determine the 
parameters {p r } for a higher-order composition [35, 36]. In Eqs.(57a) and 
(57b) it needs about r « 2 m unitary operators Sa(x) (S b (x)) to compose 
the unitary operator f2m-\{ x ) (f2m-i( x )) f° r sma U m [35]. Note that there 
is difference of the product order of the operators A and B in the opera- 
tors f2m-i{ x ) an d firn-i( x )- The difference is necessary to construct further 
the 2m— order cross interaction with the operators fztn-i{ x ) an d f 2 m~i( x )- 
Again one uses the above symmetric BCH formula to compose further the 
2m— order cross interaction: 

exp(±x(A + B) + \x 2m+l B 2m+1 + ...) exp(-x(A + B) - x 2m+1 A 2m+1 + ...) 

x exp(±x(A + B) + \x 2m+1 B 2m+l + ...) 

= exp[-x 2m+1 (A 2m+1 - B 2m+1 ) + ...]. (58) 
Obviously, here the parameter x does not take a small value so that the 
2m— order cross interaction A 2m+ i or B 2m+ i could be large enough to de- 
rive the conversion of the initial density operator to the multiple-quantum 
coherence in a spin ensemble. This is different from the conventional decom- 
position of exponential operators [33, 35, 36] where the parameter x is always 
taken as a small quantity. 

Here leave two questions of the quantum search NMR multiple-quantum 
pulse sequence in the future. One is how high order cross interaction is nec- 
essary to excite effectively the multiple-quantum coherence in an n— qubit 
spin ensemble. Another is how to extract experimentally the information of 
the marked state from the multiple-quantum spectrum excited by the higher 
order cross interaction. 



4. Discussion 

The quantum mechanical unitary dynamical description for the Turing 
machines was first given by Benioff [39] in 1980. This quantum Turing ma- 
chine model has not proved to be superior to the classical counterparts, 
as pointed out by Deutsch [2]. However, the quantum mechanical unitary 
dynamics plays an important role in computational speedup in ensemble 
quantum computation in contrast to in the pure-state quantum computa- 
tion. A direct result of the quantum mechanical unitary dynamics is that 
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the quantum search problem and the parity problem could be solved in poly- 
nomial time on an NMR ensemble quantum computer, while these problems 
are NP problems and can not be solved efficiently in classical computation. 
Therefore these works really move the frontier between the solvability and 
intractability of complex problems. There is difference between ensemble 
quantum computation based on the quantum mechanical unitary dynamics 
and the pure state version. For example, in the paper it has been shown 
that the Grover algorithm is an ineffective iterative sequence used to am- 
plify the output signal of the NMR ensemble quantum computation, while it 
can achieve a quadratic amplification for the amplitude of the marked state 
in the pure-state quantum search. The reason for this is that the iterative 
sequence of the Grover algorithm can amplify quadratically the relative am- 
plitude of the marked state and suppress any other states in the quantum 
system, but the detectable NMR signal intensity in a spin ensemble is not 
only proportional to the relative amplitude of the marked state but also the 
ensemble-averaging magnetization strength of the spin system. The ensemble 
quantum computation based on the quantum mechanical unitary dynamics 
allows the separation each other between the initial input state and the quan- 
tum circuit in a quantum algorithm. This makes it possible for the highly 
mixed states of the spin ensemble to be used in ensemble quantum computa- 
tion. This also shows that there is difference between new ensemble quantum 
computation and the classical NMR quantum computation working on the 
concept of the pseudopure state or the effective pure state. 

In the paper it has been shown that the quantum search problem could 
be solved in polynomial time on an NMR quantum computer. On the ba- 
sis of the quantum mechanical unitary dynamics the NMR multiple-quantum 
spectroscopy has been developed to solve experimentally the quantum search 
problem. In the quantum search sequence based on the NMR multiple- 
quantum spectroscopy the information of the marked state is first loaded on 
the unitary evolution propagator that excites the multiple-quantum coher- 
ence from the initial density operator, i.e., a highly mixed state, of the spin 
ensemble. This propagator is usually constructed with the oracle unitary 
operation and oracle-independent unitary operations. The NMR multiple- 
quantum spectra that carry the information of the marked state then are 
obtained experimentally. Finally one could extract the information of the 
marked state from the multiple-quantum spectra by the conventional NMR 
data processing. Although number of the multiple-quantum coherence or 
transitions increases exponentially as the qubit number in a spin ensemble, 
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number of the multiple-quantum spectral peaks can be controlled and can 
be as low as (2n + 1) in the multiple-quantum spectra for an n— qubit spin 
system. This ensures that at least some lower-order multiple-quantum peaks 
do not reduce exponentially as the qubit number as long as the initial density 
operator is transferred effectively to the multiple-quantum coherence by the 
oracle iterative sequence of the quantum search sequence. The complexity of 
the NMR multiple-quantum experiments is therefore mainly dependent on 
the oracle iterative sequence exciting the multiple-quantum coherence and 
the extraction of the information of the marked state from the multiple- 
quantum spectra. The NMR multiple-quantum spectroscopy could provide 
a more convenient and larger space to design new ensemble quantum compu- 
tation sequences or experiments. It also could provide possible approaches 
to enhancing the output NMR signal of ensemble quantum computation and 
extending the realizable size of the NMR ensemble quantum computation. In 
particular, in the frame of the multiple-quantum spectroscopy the cross in- 
teraction between oracle-dependent and oracle-independent interactions has 
been proposed to enhance the output NMR signal of the ensemble quantum 
computation. In comparison with those NMR quantum search algorithms or 
experiments [11, 17, 40] the quantum search sequence based on the quantum 
mechanical unitary dynamics and the NMR multiple-quantum spectroscopy 
would be more powerful. 
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